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FOREWORD 

This report presents results of work conducted by Boeing Computer Services 
Company under NASA Contract DEN3-42, "An Expanded System Simulation Model for 
Solar Energy Storage." The work was a continuation of the development effort 
for SIMWEST (Simulation Model for Wind Energy Storage). The current version of 
these codes also includes solar-photovoltaic energy systems modeling. This 
program was conducted under the sponsorship of the Division of Energy Storage 
Systems, DOE, under the direction of Dr. G. C. Chang, and was administered by 
the NASA-Lewis Research Center Thermal and Mechanical Storage Section with Mr. 
L. H. Gordon and Mr. R. H. Beach as Project Managers. 

This report is in two volumes. 

I. Technical Report 
II. UNIVAC Operation Manual Revisions 

The Boeing principal investigator for this program was Dr. A. W. Warren. Co- 
investigators were Dr. Y. K. Chan and Dr. M. H. Dwarakanath. 
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1.0 INTRODUCTION 


The need for computer models of energy storage for solar energy systems (wind, 
thermal, heating and cooling, and photovoltaic) is generally recognized. Al- 
though a number of computer programs have been developed for specific solar 
applications, these programs generally have lacked flexibility or compre- 
hensive modeling capability. A general purpose, user oriented, and highly 
adaptable software program was needed for solar systems research and concept 
evaluation. In November of 1976, NASA-Lewis as project manager for the Energy 
Research and Development Administration, awarded a contract to Boeing Computer 
Services (BCS), to develop such a model. It is called SIMWEST (Simulation 
Model for Wind Energy Storage). SIMWEST was required to have the capability of 
modeling wind energy/storage systems utilizing any combination of five type's 
of storage (pumped hydro, battery, thermal, flywheel and pneumatic). The level 
of detail specified was consistent with evaluating long term economic feasi- 
bility and performance of wind energy systems. 

To meet these requirements BCS developed SIMWEST based upon an existing dynamic 
simulation program (EASY). It consists of a library of system components and a 
precompiler program which allows these components to be put together in build- 
ing block form. A simulation program is then used for system evaluation. 
SIMWEST was originally developed to run on the UNIVAC 1100 series of computers. 

In July 1978, BCS was awarded a follow-on contract to expand the capabilities 
of the SIMWEST computer codes. Three areas were selected for further develop- 
ment: 

• Expansion of the SIMWEST component library, to include solar-photo- 
voltaic in addition to wind-energy generation systems, 

• Extension of program availability and support for both CDC and UNIVAC 
program versions. 
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• Specification of dynamic library components and additional analyses 
for evaluation of wind energy and storage system power transients. 

An initial evaluation of the power transient requirements indicated that to 
achieve the required degree of generality and user capability was not within 
the scope of the present contract. Subsequent program effort was then focused 
on the first two areas. The expanded SIMWEST codes enable a user to model a 
large variety of solar systems utilizing; 

• Solar-photovoltaic and/or wind power generation 

• Five types of energy storage 

• User specified system and control logic 

• Utility and backup generation 

The simulation program has proven to be efficient and versatile for performing 
parametric studies and for evaluating concept feasibility. SIMWEST is now 
available and documented for use on the UNIVAC 1100, CDC 6000 and CDC Cyber 
series computers. 

This report provides an overview of the SIMWEST program and a technical discus- 
sion of the work accomplished under the SIMWEST follow-on contract DEN3-42. 
Section 1 presents a general overview of SIMWEST modeling and simulation capa- 
bility. Section 2 describes the solar-photovoltaic components and presents 
several test case models illustrating their usage. Section 3 provides test 
case results which illustrate the capability of the SIMWEST library to model 
complex wind or solar systems including storage, loads, and system logic. 
Section 4 describes the technical approach recommended for power transient 
analyses and presents details of the component models formulated for the dy- 
namic library. 

1.1 SIMWEST OVERVIEW 

SIMWEST consists of two basic programs, and a library of generation, storage, 
environmental, and load components. The first program, the Model Generation 
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Program, is a precompiler which generates computer models (in FORTRAN) of 
complex energy generation/storage systems, from user specifications using SIM- 
WEST library components. The second program utilizes the resulting computer 
model to perform cost and power utilization analysis. It handles input, 
output, integration of system dynamics, and iterates to obtain convergence of 
implicit variables. The combination of these two programs provides a powerful 
tool for analyzing alternate generation and storage system designs. 

Figure 1.1-1 shows the general organization of the SIMWEST program. In addi- 
tion to the two programs described above, there is a third which performs file 
maintenance. It is used to incorporate user supplied data for new subsystem 
models. Although the program is shown as a number of subprograms, it can be 
executed as a single batch program by supplying the model description cards and 
the control cards describing the desired analysis to be performed and the 
desired tabular and/or plotted output. 

The SIMWEST model generation and simulation programs have a number of user 
oriented features which greatly enhance the value of the codes. Some of the 
more prominent features are shown in Table 1.1-1. These features and the 
supplemental components described in 1.2 enable the user to quickly build, 
debug, simulate and interpret alternative system designs. 

1.2 SIMWEST LIBRARY 

The SIMWEST library is listed in Table 1.2-1. It is made up of six types of 
components: environmental, generation, load, logical, storage and supple- 
mental. The two character mnemonic names are used to identify components in 
the. users model. 

The degree of detail in the component models is based upon two design criteria. 
First, all models should contain sufficient detail to simulate all physical 
characteristics and constraints having significant impact on system cost ef- 
fectiveness. Second, the models should be designed to minimize computer time 
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Figure 1,1-1 SIMWEST Program Organization 







Table 1.1-1 SIMWEST User Oriented Features 
MODEL GENERATION PROGRAM 

• Simplified Component Connections 

• Availability of all Input Parameters for Connection 

• Fortran Insertion Capability Between Components 

• Line Printer Schematic of User's Model Provided 

• Automated Naming of Parameters and Variables 

• Built-in Diagnostic Capabilities 
SIMULATION PROGRAM 

• Free Field Data Inputs, Including Tables 

• Diagnostics on Data Inputs 

• Default Values Assigned to Unspecified Parameters 

• Optional Levels of Line Printer and Diagnostic Output 

• Multiple, Back-to-back Simulation Capability 

• Printer Plotter Output of Time Histories and Crossplots 
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Table 1.2-1 SIMWEST Library Components 


ENVIRONMENTAL 


BATTERY STORAGE 


HIND 

WD 

INVERTER 

IV 

AMBIENT TEMP 

TP 

RECTIFIER 

RE 

TMY WEATHER TAPE 

ED 

BATTERY 

BA 


\ 

ADMITTANCE 

AD 

WIND POWER GENERATION 




TURBINE/GENERATOR 

WP 

FLYWHEEL STORAGE 


WIND TURBINE 

HT 

AC MOTOR 

MO 

FIXED RATIO TRANSMISSION 

GR 

VARIABLE RATIO TRANSMISSION 

TR 

AC GENERATOR 

GE 

FLYWHEEL/CLUTCH 

FL 

SOLAR POWER GENERATION 


HYDRO STORAGE 


SOLAR ORIENTATION (TRACKING) 

SO 

HYDRO: PUMP 

PU 

FLAT PLATE COLLECTOR 

FP 

HYDRO TURBINE 

HT 

FOCUSING LENS COLLECTOR 

FO 

HYDRO STORAGE 

HS 

PHOTOVOLTAIC ARRAY 

PV 

PNEUMATIC STORAGE 


UTILITY GENERATION 


COMPRESSOR 

CO 

UTILITY 

ur 

TURBINE 

TU 



ADIABATIC HEAT EXCHANGER 

HX.HY 

LOGIC 


BURNER 

BN 

POWER DIVIDER 

PD 

PNEUMATIC STORAGE 

cs 

POWER ACCUMULATOR 

PA 

THERMAL STORAGE 


PRIORITY INTERRUPT 

PI 

STORAGE VESSEL 

TS 

SWITCHES SW.SX.SY/SZ 





IM 


SUPPLEMENTAL 


ELECTRICAL LOAD 

LO 

SATURATION 

SA 

THERMAL LOAD 

TL 

RANDOM NUMBER GENERATOR 

RN 

V i ' 


TEST FUNCTIONS 

af 

SUPPLEMENTAL 


TABLE LOOKUPS 

FU.FV 

ARITHMETIC ELEMENTS 

MA.MBjMC 

TRANSFER FUNCTIONS IT,UUL,TF 

COST MONITOR 

CM 



HISTOGRAM 

HG 



TAPE READ 

TA 



TIME CONVERSION 

TI 
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and required user specification. It is assumed that a Sl'MWEST simulation might 
cover a time span of one year. Thus, from a computer run time and economic 
impact point of view a simulation step size of between 15 minutes and one hour 
was established as a design goal. 


As a result of the above design criteria, many physical components, such as the 
electrical components, were modeled mainly in terms of power flow and steady 
state response. This level of detail is consistent with a 15 minute time step 
and with the concept that important transients are on the time scale of demand 
curves or weather patterns, i.e. , an hour "or more, rather than on the time 
scale of electric motor transients of a few seconds. If short time transients 


were to be modeled, additional detail would be required in the component models 
which would greatly increase the user's task of specifying the model. Further, 
the simulation time step would have to be reduced and computer runs would be 


much costlier. 


The environmental components listed in Table 1,24 simulate environmental con- 
ditions. In the present SIMWEST library a user can generate wind speed and 
ambient temperatures, or can use selected Inputs from the recorded weather and 
insolation data on the Typical Meterological Year (TMY) tapes for one of 26 
U,$, locations. These variables are generally used as inputs to physical 
components. 

The generation components consist of wind generation, solar-photovoltaic and 
utility routines, The wind turbine-generation components are fairly simple 
models for computing the power output of a conventional, horizontal axis wind 
machine given basic machine parameters. The solar-photovoltaic components are 
somewhat more sophisticated, especially in the collector thermal analysis, and 
have a number of modeling options which a user may employ, e.g. , active or 
passive cooling. 

The storage components encompass such things as motors, generators, trans- 
missions, and flywheels. These components model actual physical hardware 
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which might be used in a wind or solar energy system. The selection of the 
particular SIMWEST library set of storage components was based on the require- 
ment that it be capable of modeling the five types of energy storage systems 
mentioned previously: thermal, flywheel, battery, pumped hydro and pneumatic. 

The load components in the SIMWEST library are used to simulate various types 
of power demand. They also monitor how well the system meets the simulated 
demand and compute the value of the energy delivered to the load. Like the 
environmental components, these components may be computed from actual mea- 
surement data or from randomly generated data based on user furnished load 
profiles. 

The library's logical components are the power dividers, power accumulators, 
switches and priority interrupts. Although physical hardware or logic devices 
could be built to serve the function of the logical components, they are not 
meant to represent any particular existing hardware. Instead, they are ideal- 
ized components that allow the user flexibility in modeling a wide variety of 
system and control logic for operational evaluation of energy storage systems. 
In practice, the control function might be performed by a control room operator 
using a predefined control strategy or by use of a process computer. 

Finally, the supplemental components include such things as the tape read, the 
histogram and the cost monitor. These components serve to help the user run 
the simulation and analyze its results. 

1.2.1 Storage Subsystems : : 

Figures 1.2-1 through 1.2-5 give possible configurations of the five types of 
storage subsystems which can be modeled with the present SIMWEST library. For 
illustrative purposes the number of variables shown passed between components 
is limited. A description of the variables being passed is given in Table 
1 . 2 - 2 . 
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Figure 1.2-1 Pneumatic Storage Subsystem 
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Figure 1.2-3 Flywheel Storage 
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Figure 1.2-4 Battery Storage 
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Figure 1.2-5 Thermal Storage 
Table 1.2-2 Partial List of Component Inputs and Outputs 



SYMBOLS 

1 

pj 

POWER 

RE 

POWER REQUEST 

MP 

MAXIMUM POWER 

RS 

ROTOR SPEED 

T 

TEMPERATURE 

TA 

AMBIENT TEMPERATURE 

M 

MASS FLOW RATE 

H 

RESERVOIR HEIGHT 

LD 

THERMAL LOAD DELIVERED 

wv S' 

WIND VELOCITY 

GR 

GEAR RATIO 

EF 

EFFICIENCY 

INT 

INTERRUPT FLAG 

PR 

PRESSURE 

PS 

PRIORITY SEQUENCE 

WY 1 

i i . . V ■ 

WEEK OF YEAR 

DW 

DAY OF WEEK 

TO 

TIME OF DAY 

SP 

SURPLUS POWER 
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A total energy system will generally be made up of elements from a number of 
different subsystems (see Figure 1.2-6). In addition, the SIMWEST program can 
be used for models which include networks of storage subsystems of the same 
type or a network of wind or solar generators. 

1.2.2 Logic Components 

The capability for modeling complex system control logic is provided by the 
power divider, power accumulator and priority interrupt components. Both the 
divider and accumulator operate on a priority basis. The priority interrupt is 
used by other system components to change the priority setting of the divider 
and accumulator. 

The power divider has one input power port and four output power ports (not all 
output ports need be used for a given simulation). The divider also has an 
input request associated with each of its output ports. A power request 
originates with a component which is directly or indirectly connected to an 
output port. The user specifies priorities of either 0, 1, 2, 3, or 4 to be 
associated with each of the output ports. If the input power exceeds that 
requested of the port with highest priority (priority 1) then the excess power 
goes to the port with the next priority. This process continues until either 
all power is distributed or all requests of non-zero priority ports are met. A 
port with zero (0) priority does not receive power. Such ports are included to 
model backup or switch operated components. In these situations, the connected 
component would change the zero priority setting of the power divider by use of 
a priority interrupt. Two or more ports may be assigned the same priority in 
which case the user may specify weights to be associated with each port. Then 
if there is not enough power available to satisfy all requests of equal prior- 
ity, the power is divided between them in proportion to the user specified 
weights. 

The power accumulator is similar to the divider except that instead of distri- 
buting power from a single input port among four output ports, it accumulates 
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Figure 1.2-6 Example of Power Divider and Accumulator Use 
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power from four input ports and sends it out through a single output port. The 
power accumulator accepts power requests from the downstream component and 
al Ideates requests to each of its input ports in order to service the down- 
stream component. • 

An example illustration of the use of power dividers and power accumulators is 
given in Figure 1.2-6. It is seen that power from the turbine/generator is 
distributed with highest priority (priority 1) going to the power accumulator 
that services load 1. Since the power accumulator servicing load 1 has its 
priority 1 input port connected to the power divide^, it will try first to 
satisfy load 1 from the turbine/ generator and then from the utility. If the 
power divider satisfies load 1 and there is power left over, it will be used to 
satisfy the request from the battery. Finally, if the battery is full or if 
its charging rate is met, then the excess power goes to the flywheel. The 
battery also has a priority zero connection to the utility. If the battery 
remains in a discharge state for more than a specified amount of time, it can 
change the utility priority (from 0 to 1) to receive needed power. 

Also in Figure 1.2-6, we see that load 2 prefers to draw power from the 
flywheel before turning to the battery. This configuration tends to keep the 
flywheel as discharged as possible, using it primarily as a means to absorb 
large influxes of power. 

1.3 SIMWEST OUTPUT 

There are three basic forms of SIMWEST output to facilitate the analysis of 
wind and solar energy storage systems; line printer: plots, histograms of system 
variables and time sequenced output of variable values. Each SIMWEST library 
component is associated with a number of output variables. Prior to simulating 
a given system the user may select any of these outputs for plotting or tabular 
output. For example, he may want to plot the energy of pneumatic storage as a 
function of time and/or as a function of temperature. If the user wants a time 
sequenced listing of all variable values, he may specify the time step between 
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printouts. The listing of all variables has proven to be a useful tool in 
understanding the performance of the system under consideration and a valuable 
aid in validating the system design. 

■m 

SIMWEST also provides a special output which computes life cycle and levelized 
energy costs per kwh. This output is produced by the cost monitor component 
and is illustrated in Figure 1.3-1. The levelized energy costs are based on 
energy delivered to the loads during a simulation and forecasted! to a full 
years' system operation. This output permits direct comparison of capital and 
energy costs for alternative system configurations, enabling a user to perform 
economic trade studies and system sizing. 

1.4 TESTING 

Reference [ 1 ] describes two simulation studies which were used to test the 
original SIMWEST program. Section 3 describes the NASA-Lewis approved simula- 
tion studies for the expanded SIMWEST program. These studies provide an 
excellent test and illustration of the program's capability to model complex 
wind/solar energy systems. 

Prior to performing the simulation studies and throughout its development the 
SIMWEST program was systematical ly tested. First components were grouped into 
simple systems and simulations were performed. During these simulations sys- 
tem parameters were driven so as to force the individual components through 
every normal program path and to assure that all component outputs assume a 
wide range of values. The number of components and the number of ways they can 
be connected makes it impossible to exercise every combination. However, the 
subsystem groupings that were used were representative of the expected program 
usage. Section 2 describes some of the test cases for the solar-photovoltaic 
generation components. 

The test cases and simulation studies revealed that the code is reasonably 
efficient for system parameter studies. Even on very complex systems, such as 
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SOLAR/WIND ENERGY STORAGE COST SUHMARY 


20 YEAR LIFE CYCLE 


i YEARLY SYSTEM COSTS 

I 



CAPITAL COST 

(INCLUDING FIXEO CHARGES) 

526* S 


FIXED 0 ♦ M COST 

107* S 


OPERATING ♦ FUEL COST 

14* * 


TOTAL 

*48. S 

ENERGY 

DELIVERED 



ENERGY DELIVERED 

7445* !C UH 

* 

• 

ENERGY COST PER KWH 

* 

86.8 MILLS * 

: • 


VALUE OF ENERGY OELIVERED 
(VALUE OF FUEL SAVED) 

372 S 


ENERGY VALUE PER KVH 

50*0 MILLS 


COST PER VALUE DELIVERED 

1.74 

LOAO FACTOR 



PERCENT OF LOAD SUPPLIED 
BY total solar SYSTEM 

100 *0 


PERCENT OF LOAO SUPPLIED 
BY UTILITY 

0*0 


PERCENT OF SOLAR ENERGY 
SURPLUSEO 

0*0 


COST TO MEET LOAD 
(SOLAR ♦UTILITY) 

86*8 MILLS 


Figure 1.3-1 Cost Monitor Output For Fresnel Lens Model 

(Example Only) 
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represented by the NASA-Lewis test case, convergence of logic variables was 
quite rapid. Convergence generally took place in less than six iterations per 
simulation time step. As an example, the year simulations used in the NASA 
defined parameter study of Reference [1] took less than 420 CPU seconds on the 
CDC 6600. For comparison, the CPU time on the UNIVAC 1100/40 is approximately 
two to three times as great as that on the 6600, and CPU time on the Cyber 175 
is a factor of two to three times smaller than that of the 6600. 

1.5 PROGRAM USAGE 

During the testing it became clear that while the user need not be a SIMWEST 
expert or software specialist to make efficient use of the program, he should 
thoroughly think through and be familiar with the characteristics of the system 
he wants to simulate. Component models, if not carefully specified, may 
perform in unexpected ways. If the systems logic is not well thought out, the 
resulting system may be significantly out of balance and subsystems may not be 
fully utilized. The test case described in Section 3 illustrates the process 
of sizing and logic adjustment to satisfy system performance objectives. 

A number of useful procedures were developed during the simulation studies. 
First it was found that when simulating a complex system, it is best to 
separately develop and test subsystem portions of the model. This allows 
problems or unexpected results to be isolated and understood prior to the 
introduction of the more complex characteristics associated with the total 
system. 

It was found during the simulations that the use of Fortran statements in the 
model definition is very useful for creating special input to system components 
and for defining special outputs to be plotted or statistics to be printed. 
For example, Fortran statements enable the user to generate and interpret trade 
study data by computing component input parameters from user specified system 
parameters. The use of Fortran statements is simple and should be encouraged 
early in SIMWEST applications. 
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Computer simulation costs may be minimized by appropriate tradeoffs between 
run time and simulation accuracy. Run time is most directly affected by the 
integration step size, the total simulation length, and the average number of 
iterations through the model at each time step. For long duration runs, an 
hour step size is usually acceptable. Models having smaller time constants 
than the step size may be approximated by implicit steady state conditions and 
solved by iteration through the model. If a model requires many iterations for 
convergence then it may be useful to isolate the source of instability in order 
to modify or simplify that portion of the system model. It has been generally 
found in the simulation studies that use of a few seasonal weekly simulations 
is adequate to predict long term performance for system trade studies and 
design optimization. Based on the results of Section 3, four to six week long 
simulations are recommended for this purpose. 

When making a year simulation run, it is best to break it into twelve monthly 
simulations. Thus, measures of performance such as plots, histograms and 
performance statistics are available on a monthly basis. In addition to giving 
better visibility of the system performance, this helps limit the job core 
size. The twelve monthly simulations can be submitted as a single run with the 
results of a given month acting as initial conditions for the next month. The 
user only needs to submit new data cards for data which changes from one month 
to the next. 
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2.0 SOLAR PHOTOVOLTAIC COMPONENT MODELS 


The solar photovoltaic component models added to the SIMWEST library as a part 
of contract DEN3-42 are briefly described and test case results illustrating 
their use are summarized in this section. 

Table 2.0-1 summarizes the characteristics of the solar-photovoltaic compo- 
nents. The environmental data component is designed to read Typical Metero- 
logical Year (TMY) data tapes containing hourly insolation and weather data at 
26 U.S. locations. This component can also be used to read other hourly data 
tapes such as the SOLMET tapes by inputing a user specified format to the model 
generation program. The solar orientation or tracking component computes the 
sum of direct beam and global insolation on a flat plate array for fixed 
orientation and four different beam tracking options. The flat plate and 
focusing lens collector components provide detailed thermal analyses for de- 
termining average solar cell temperature. The collector models, and that of 
the solar array are based on similar models developed at Sandia Laboratories 
for the SOLCEL program. (Reference [4] ) . The array component model is a sim- 
plified model based on scaling the characteristics of a single solar cell. 
Array voltage can either be user specified or determined by a maximum power 
tracker. It should be observed that the above components are coded in SI 
(metric) units, whereas most of the SIMWEST components are coded in English 
units. This is generally not a problem since there are at most only a few 
interconnection variables between the solar-photovoltaic generation components 
and other SIMWEST components, and these are easily converted using arithmetic 
components. 

The TMY data tapes are currently the best environmental data sources available 
for simulating typical yearly solar energy system performance. These tapes 
were extracted from SOLMET data tapes containing rehabilitated hourly solar 
and meterological observation data over a period of many years at each observa- 
tion site. Each Typical Meterological Year was created by statistical selec- 
tion of a typical meterological month for each calendar month in the long term 
data base and catenating the 12 months to form a TMY. All of the TMY data files 
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TABLE 2,0-1 SOLAR-PHOTOVOLTAIC COMPONENTS 


COMPONENTS SYMBOL 

• -ENVIRONMENTAL ED 

• SOLAR ORIENTATION SO 

• FLAT PLATE FP 

COLLECTOR T 

• FOCUSING LENS FO 

COLLECTOR 

• PHOTOVOLTAIC ARRAY PV 


PURPOSE 


READ DOE SOLAR INSOLATION AND WEATHER 
DATA-TYPICAL METEOROLOGICAL YEAR TAPE 

SOLAR INSOLATION ON TILTED FU\T PLATE 
ARRAY. (FIVE OPTIONS) 

FLAT PLATE THERMAL MODEL WITH FLUID 
AND PASSIVE COOLING OPTIONS 

FRESNEL LENS THERMAL MODEL WITH FLUID 
AMD PASSIVE COOLING OPTIONS 

CONVERTS SOLAR INSOLATION TO D.C. 
ELECTRICAL POWER. MAXIMUM POWER TRACKER 
OR USER SPECIFIED VOLTAGE 


are available for use by a SIMWEST user. He thus has access to a high quality 
environmental data base for solar energy simulations and system analyses. 

2.1 PHOTOVOLTAIC MODEL TEST CASE 

The input data for the photovoltaic model test case is shown in Figure 2.1-1. 
The purpose of this model is to obtain characteristic current voltage curves 
for the default solar array parameters. Fortran statements are used in the 
model generation data to let the terminal voltage range between 0 and 204 volts 
for solar insolation values of 5, 20, and 50 suns (1 sun * 1000 w/m 2 ). Cell 
temperature is specified at 25°C for the first simulation and 55°C for the 
second. Figure 2.1-2 shows the current voltage curves and Figure 2.1-3 shows 
power voltage cross plots at the lower cell temperature and for the three solar 
insolation levels. These curves verify the physical characteristics of the 
solar cell model. It may be noted in these figures that current and output 
power become negative when the specified voltage exceeds the array open circuit 
voltage. Individual cell characteristics may be obtained by dividing voltage 
by 300 (default number of cells in series) and by dividing current by 500 
(default number of cells in parallel). 

2.2 FLAT PLATE COLLECTOR MODEL 

The input data for the flat plate model test case is shown in Figure 2.2-1. The 
purpose of this model is to illustrate water and wind cooling of the collector 
and to test the tracking options of the orientation component SO. There are 
six 1-1/2 day simulation runs. The first run uses water cooling (CM0FP=2), a 
single glass cover over the front plate and insolation on the back. The second 
run uses passive cooling (CMQFP=Q), no plate insolation and fins on the back to 
cool the collector. In the first two runs, the collector is tilted and has a 
fixed, southward facing orientation (MO $0=1). The last four runs are similar 
to run 2 except different tracking options are utilized. 

The model schematic produced by the model generation program is shown in 
Figure 2.2-2. The component TI is used to furnish time of day and day of year 
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MODEL DESCRIPTION PHOTO-VOL TA IC 
LOCATION=l 1 TI 

k— FORTRAN— STATEMENTS 

ST PVsSOOO 

JFCDY TI,GT,U5)ST PVs20000 

IHDY-TI-, GT t 2t5 ) S T-PVsSOOOO 

VT PV=8,5*TD TI 
LOCATXONB53 PV 

— EflD-OF— MOOEt 

PRINT 


PHOTO-VOLTAIC CURRENT VOLTAGE CURVES 


a) Model Generation Input Data 


PARAMETER VALUES 

cycles=o,to ti=o 

— DLI M E S = 5 0 

TC PV=2b 
RC PV=1 

—PRIM TER— PLOTS# DISPLAY 1 
V PV # VS# TIME 
I PV#VS#V PV 


P RV#VS#TIME 

TINCb,5,TMAXs72#PRATEs 24#PRINT CONTROLS# INT M0DE=3,0UTRATEst 

—T-HLEs PHOTO* VOL TA I C— CE LL— C U R R EN T— V 0 IrT A G E~~C U R VLS 

SIMULATE 


PARAMETER VALUES 
— 4C— P-V«55 

simulate 



b) Simulation Program Input Data 



Figure 2.1-1 PV Test Case Input Data 
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Figure 2.1-3 Solar Array Output Power Versus Voltage 










'model description FLAT PLATE test case 

LOCATIONS U TI 

ikOCA-TlDN = 35 ED INPUTS«TI 

IL0CATIONS53 SO INPUTSsTI #ED(X1=SB#X2=8T) 

ffLOCATION=57 FP INPUTS=SO,EDCX«=WD,X3 sTA) 

■ fiND-OF-fiODEL 

PRINT 


a) Model Generation Input Data 


parameter values 

CYCLE S=2, 01 f TO T I a36 # TF IFPs 1 0 # TFOFPs30# MFMP Ps , 02, CMOF P«2 # NG FP=I# 

OL-INESsSO t 

Hi! F P s f 0 t 

XW|Fpsi,CL FP=2,NT FPS10/CC FPM000,CM FPs l 0,CPOFP= , 0 i #LA $Os29,733# 

■TL— S0=29 t -733, AA-SOB2 <■ 

PRINTER PLOTS# DISPLAY! 

TLTSO# VS# TIME 

4c— FP , V 3 #-T-I M E 

X2 ED, VS# TIME 
PI FP# VS# TIME 

-T'INC = v5#-TMAXs36,PRATEs6#PRINT— CONTROL-s3#INT-MODEs3#OUTRATEsi 

TITLEsPLAT PLATE COLLECTOR TEST CASE 
SIMULATE 

-PA R A M F.-TE R -V A LUES 

CMOFPsO,HI FP=i,E9#FIRFPs« 

SIMULATE 

-P A R A h E-T E R -V ALUE S — 

MO SUs2 

simulate 

-PARAMETER— VALUES 

MQ SQsJ 
; SIMULATE 

f -PARAMETER— VALUE 8 — 

MO SO=4 
f i SIMULATE 

-FAR-AME-TE-K-VAL-UES : 

MO SQs5 
; • SIMULATE 

I •• ~ - ■ . " ■- . ' , ■ .. ■■ 


b) Simulation Program Input Data 
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Figure 2.2-1 Flat Plate Collector Model Input Data 
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information to SO and to the TMY read component ED. ItD supplies direct beam and 
global insolation to SO, and ambient temperature and wind speed to the col- 
lector component FP. Based on collector orientation, SO supplies solar inso- 
lation incident to the array, collector tilt angle, and tracking power to FP. 

Typical results of the flat plate model runs are shown in Figures 2.2*3 through 
2.2-5. Figure 2.2-3 shows the global horizontal insolation obtained from ED 
during the 36 hour simulation period. The data period was during mid-winter 
and the daily peak levels are thus low to moderate. The array tilt angle daily 
pattern for horizontal E-W axis tracking is shown in Figure 2.2-4. At noon the 
array is oriented normal to the sun's incident rays and thus maximizes the 
insolation gathered during the mid-day peak. The tilt angle approaches 90° as 
the sun approaches the horizon, and remains fixed at 90° overnight. Comparison 
of the solar insolation peaks with the various tracking options showed that 
horizontal E-W axis tracking gave the best results of the single axis tracking 
systems, and was only slightly inferior to two-axis beam tracking. Solar cell 
temperature for this case is shown in Figure 2.2-5. The cell temperature is 
within a few degrees of ambient most of the day and rises in mid-day propor- 
tional to the solar insolation received. The results with water cooling are 
quite similar. 

2.3 FRESNEL LENS COLLECTOR MODEL AND INCREMENTAL COSTS 

The input data for the Fresnel Lens test case is shown in Figure 2.3-1. The 
purpose of this model is to illustrate a Fresnel Lens collector model with 
thermal fluid loops for collector cooling and for solar heating. Three week 
long simulations are used to demonstrate incremental cost calculations for 
subsystem economic design. A variable speed pump is assumed for the collector 
fluid loop with the flow rate adjusted so that the outlet temperature is 5°C 
greater than the inlet. The collector consists of a rectangular grid of 120 
Fresnel lenses each of which focuses solar radiation on a 5 x 5 array of solar 
cells. Excess thermal energy is conducted to a heat sink surface and then 
dissipated by natural convection, radiation and heat exchange to the coolant 
fluid. The collector parameters are chosen for a lens concentration ratio of 
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Figure 2„2-3 Global Horizontal Radiation Versus Time 
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Figure 2.2-5 Solar Cell Temperature Versus Time 
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MODEL DESCRIPTION 

FRESNEL LENS COLLECTOR WITH THERMAL STORAGE AND LOAC 

LOCATI ON=ll 

TI 


LOCAT I ON- 71 

ED 

INPUTS=TI 

LOCAT IjON=45 

MA 

INPUTS=TS(T=FIN) 

FORTRAN STATEMENTS 


TFOFO 

= FO MA ♦ 

5 . 

LOCAT ION-33 

FO 

INPUTS-ED(X1=ST, X3=TA,X4=WD) ,MA{FO=TFI) ■ M 

LOCATI 0N=73 

PV 

INPUTS=ED(X1=ST) ,FO 

LOCATIONS 47 

TS 

INPUTS=FOCP,l=P),TL 

LOCA T I ON=27 

TL 

INPUTS=TI ,ED(X3=TA> 

LOCATI ON=77 

LO 

INPUTS=PVCP=P,l,P=LO,l) 

LOCAT I ON=79 

CM 


END OF MODEL 



PRINT 



a) Model 

Generation Input Data 


TITLE=FRESNEL LENS COLLECTOR ( INCREMENTAL COST COMPUTATION) 
i PARAMETER VALUES 

CYCLES=4.01,T0 T 1=0 ♦ CM0F0=2, CW F0=3.75,CL F0=3.9* DLINES=5Q 
fr NL F0=12Q,NT F0= 24 . MFMF 0 =0 .3 , CC F0=6.,CM FO=50,HI FO=,01tRC F0=.06 
U TS TS=5,0H TS=. 00879, PD TS=12,LE TS=30,NU TS=.01,NC TL=0.2 
i f Cll M A~* 55556 »C2 MA=-17.7778» COPFO=0*5 
; 'f cci PV=10a,CM PV=50,LE TS=30,CR CM=15,LE CM=20 
1 AA P V=0 • S , NS PV = 60dtNP P V=5, RAPPV=1.3 
[ P VE LO=.05,VE TL=« 05 
!■ f TABLE* HT TS = 4 

T • 00879, .025491,-047371, .064072 
::i 90,147,147*204 
M TABLE, TL0TL=4 
ft -10,0,10,25 
- 4, 2,1.5, 1 
it TABLE, TWTTL=4 
if 0,6,13,24 
• 4, 1,1, .4 

It PRINTER PLOTS, DISPLAY1 
it RE TL, VS, TIME 
I E TS, VS, TIME t 
1 PI FO, VS, TIME 
5 FMDFO,VS,TIME 
DISPLAY2 
TC FO, VS, TIME 
P PV , VS, TIME 
FO HA, VS, TIME 

INITIAL CONDITIONS=E TS=80 

TINC=« 5 , TMAX=168 , PRA TE=l 2* PR INT C0NTR0L=3, INT M0DE=3 , OUTR ATE = 1 
SIMULATE 

[ PARAMETER VALUES* TS TS=5.5 
! SIMULATE 

PARAMETER VALUES 

TS TS=5.,NL F 0=126, CW F0=3.94,AA PV=0.63,NS PV=630 
SIMULATE 

j b) Simulation Program Input Data 

Figure 2.3-1 Fresnel Lens Model Input Data 
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25 and series connection of the output from each array. At maximum output the 
array collects about lOkw of solar radiation and produces about 1.7kw of 
electrical power. 


The model schematic produced by the model generation program is shown in Figure 
2.3-2. The collector thermal loop is formed by the connections between the 
collector FO, the thermal storage TS, and the multiply and add component MA. 
The MA component is used to convert the thermal storage outlet temperature from 
degrees fahrenheit to degrees centigrade. The output temperature from MA is 
supplied as the inlet temperature to FO. The total thermal power gathered by 
the coolant fluid is computed in FO and supplied to TS. Similarly, the thermal 
load fluid loop is represented by a power request from the load component TL to 
TS and by thermal power delivered from TS to TL. The electrical output of the 
array is computed by PV and supplied to a load component LO which monitors the 
electrical energy collected. 


Results of the. first week simulation run are summarized in Figures 2.3-3 
through 2.3-6. The weather was fairly constant during this run and solar 
insolation was fairly strong all week. Figure 2.3-3 shows that with water 
cooling cell temperature was held to less than 70°C at peak insolation. In 
fact, about 60 % of the solar energy incident on the array is exchanged to the 
coolant fluid during peak insolation. The electrical output of the array is 
shown in Figure 2.3-4. The fluid flow rate of the pump and thermal energy 
collected exhibit very similar daily patterns. The thermal load for this week 
is shown in Figure 2.3-5. This load is dependent on both time of day and 
ambient temperature which yields the complex load pattern shown. Figure 2.3-6 
shows the temperature of the thermal storage vessel resulting from the col- 
lector and load thermal loops. The daily cycles are predominant with the 
periods of strong insolation providing sufficient energy to satisfy the load 
and compensate for thermal losses. Average load is fairly well matched to 
solar generation during the week since the temperature remains within a 15° 
channel and does not have an apparent trend away from this range. 
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Figure 2.3-2 Fresnel tens Model Schematic 
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79/03/09. 12.43.33. SIMULATION DISPLAY 2 CASE NO. 1 

Figure 2.3-3 Solar Cell Temperature for One Week Simulation 
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Figure 2.3-4 Photovoltaic Array Output for One Week Simulation 
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Figure 2.3-5 Thermal Load Demand for One Week Simulation 
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Figure 2.3-6 Thermal Storage Temperature for One Week Simulation 


One of the most important measures of performance for a solar energy system is 
the level ized cost of energy, i.e., the ; life cycle cost to produce one unit of 
usable energy including generation, storage, transmission and conversion sub- 
systems. Energy cost may be used to size components and select most promising 
system alternatives, i.e., minimum energy cost is used as a selection or 
optimization principle. Although SIMWEST does not provide user optimization 
capability, optimal sizing of a few key parameters, such as the ratio of solar 
to utility generation and the size of storage relative to generation, is 
possible and may be accomplished quickly using the concept of incremental 
energy cost. The idea is to compute the incremental change in levelized energy 
cost per incremental change in capital cost, for the system parameters of 
interest. Given an initial system configuration and M sizing parameters to be 
selected, optimization proceeds as follows: 

1) Perform M+l back to back simulations to compute the cost and energy 

performance of the baseline configuration and M incremental configur- 
ations from the baseline. 

2) Calculate the incremental energy costs for each parameter variation. 

Then select a new baseline configuration. Since the incremental costs are 
equal at the minimum cost point, increase or decrease the sizing param- 
eters so as to equalize the new baseline incremental costs. 

3) Go to 1) and continue adjusting subsystem parameters until either a per- 

formance limit is reached or until the incremental costs of the remaining 
parameters are equalized. (If two incremental costs are unequal, one can 
always lower the system energy cost by increasing the subsystem with the 
smallest incremental cost at the expense of the other subsystem.) 

This procedure is recommended as more efficient and economical than using a 
series of parametric trade studies for subsystem optimization. 
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The process of computing incremental costs is illustrated for the Fresnel Lens 
model. In the first simulation the baseline system performance and costs are 
computed. The second simulation differs from the first in th-xt thermal storage 
capacity has been increased by 10%, and the third simulation differs from the 
first in that the solar collector and photovoltaic array area have been in- 
creased by 5%. Table 2.3-1 summarizes the incremental cost and simulation 
results for these runs. Column 1 shows the initial capital cost of the 
baseline system and the incremental 1 capital costs for the thermal storage and 
solar array increases. (These costs are meant to be illustrative rather than 
representative.) Column 2 shows the results of a 20 year levelized cost 
analysis of the three systems, including maintenance and operating costs, 
e.g., the change in thermal storage increases costs by $9.10 per year. Column 
3 shows the energy delivered to the loads in a year as estimated from the one 
week simulations. (Note: the change in storage capacity lowers the average 
coolant temperature, thus increasing output power.) Column 4 shows the level- 
ized energy costs of the baseline system and of the increments in storage and 
generation. This column shows that the levelized energy cost will decrease as 
thermal storage or generation are increased, and that thermal storage is under- 
sized relative to generation since a fixed $ increase in storage will lower the 
system energy cost more than the same $ increase in array area. Column 5 shows 
the % change in levelized energy cost given a 1 % increase in capital invest- 
ment. This column contains the same basic information as column 4 but provides 
a better quantitative measure of the economic value of increased storage capac- 
ity. 
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Table 2.3-1 Incremental Cost Calculations 



CC 

LC 

ED 

EC 

Baseline 

7392. 

-j i 

1272. 

7829. 

16.2 

10% Inc. in Thermal 

61. 

9.10 ; 

110.5 

8.2 

5% Inc. in Solar 

319. 

47.90 | 

365.0 

13.1 


NOMENCLATURE: CC = Initial Capital Cost in $ 

LC = Levelized Total Cost/Yr. in $ 

= Capital Cost*Life Cycle*Charge Rate + 
Maintenance Cost + Operating Cost 

ED = Useful Energy Delivered/Yr. in KWH 
= Electrical Load + Thermal Load + 

Net Change in Thermal Storage 

i 

EC = Levelized Energy Cost in <£/KWH 
= LC*100/ED 

NIC = Normalized Incremental Costs 

= % Change in EC Per % Change in CC 
= (ALC/LC - AED/ED)/(ACC/CC) 
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3.0 HYBRID SOLAR TEST CASE 

This section illustrates the capability of the expanded SIMWEST library for 
modeling complex wind or solar systems including storage, loads, and system 
logic. A test case model was chosen utilizing: 

• Typical Meteorological Year (TMY) environmental inputs 

• Hybrid generation - wind turbine plus solar array 

• Battery, flywheel and heat storage systems 

• Three user specified loads plus control loads 

• Complex system logic including utility tie-ins. 

The model and data inputs were debugged and roughly sized using week long simu- 
lations. Four seasonal week long simulations were used to balance generation 
to loads and power distribution to storage. Full year long simulations were 
then run to measure typical year system performance. This task utilized about 
three manweeks from model inception to completed year long simulations. 

3.1 SYSTEM DESCRIPTION 

An overview schematic of the hybrid solar test case model is shown in Figure 
3.1-1. Each box or function is represented by one or more components from the 
SIMWEST library, and the arrows between boxes show information and power flow 
through the model. In this model , all the power generated by the wind turbine 
and solar array are distributed to the storage devices with any remaining power 
surplused to a utility. The DC load is met from battery storage, the control 
loads from flywheel storage, and the heating load from thermal storage. The AC 
load has been segmented into baseload and peak load portions to illustrate 
utility type unit commitment to load demand. The baseload is met from battery 
storage and the peak load from flywheel storage with a utility backup if either 
storage unit has insufficient capacity. Utility backups are also provided for 
the storage devices. 
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Figure 3.1-1 Hybrid Solar Test Case Model 



































The boxes with numbers in Figure 3.1-1 represent power divider and power 
accumulator logic controllers. The flow of power in these multi-port com- 
ponents is regulated by the priority numbers with 1 = highest priority, 2 = 
second priority, 3 = third priority, 0 = no priority (used for back-up or 
switch regulated components). For example, the power divider component pre- 
ceding the battery allocates the power generated by the wind turbine and the 
solar array to the battery first (up to its maximum charging level), to the 
flywheel second, and to heat storage third, with any remaining power being 
surplused to a utility. These logic components have proven to be quite versa- 
tile for specifying and modifying model control logic. 

Once an overall system model has been defined, the user may, construct a Fortran 
system model using the SIMWEST model generation program. Figure 3.1-2 shows 
the model generation input for the hybrid model. This program produces a line 
printer schematic of all the component model interconnections which enables 
the user to quickly debug and verify his system model . Figure 3.1-3 shows the 
line printer schematic for the generation portion of the hybrid model. The 
symbols within the boxes are two or three character mnemonics for library 
components, i.e., WT represents a wind turbine component, GR a fixed gear ratio 
transmission, and GE1 the first generator component in the model (there are 
several). The four to six character symbols between boxes are the formal names 
of interconnection variables, i.e., P2 WT is the power output from the wind 
turbine into the gear box. 

The most important system characteristics of the hybrid test model are sum- 
marized in Table 3.1-1. Dodge City, Kansas TMY inputs were chosen since this 
location has both a favorable wind regime and good solar insolation. The total 
rated output of the solar array and wind turbine is 50kw, which produces an 
average output of 14kw at this location. The storage units were sized with 
sufficient capacity to accommodate most of the variability in output power due 
to changing weather fronts. Thus the battery and thermal storage units have 
about 5 days of storage, whereas the flywheel has a 1.4 day storage capacity to 
accommodate daily load fluctuations. Figure 3.1-4 shows the daily load 
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Figure 3.1-3 Line Printer Model Schematic (Page 0 Only) 
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Table 3.1-1 Hybrid Model 


G ENERATION 

• Environmental 

Dodge City, Kansas Inputs 
Mean Wind Speed = 5.2 m/s 


• Solar P-V System 

Flat Plate Array 
Horizontal EW Axis Tracking 
Passive Cooling with Fins 
Maximum Power Tracker 
Rated Output = 30 kw 
(189 m 2 Array Area) 

• Wind Turbine System 

- Horizontal Axis, Fixed Speed 
Turbine 

D.C. Generator 
Rated Output = 20 kw 
(6.4m Blade Radius) 

LOADS 1 

• AC Electrical Load 

Plumbrook, Ohio Profile 
Mean Load = 5 kw 
Base Load Power = 4.5 kw 

• DC Electrical Load 

Mean Load = 1 kw 
Peak Load = 4 kw 


• Thermal Load 

Hot Water Heating 
Mean Load = 4 kw 
Peak Load = 12 kw 


1 See Figure 3.1-3 for Load Profiles 
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System Description 


STORAGE 

• BATTERY 

Lead-Acid Battery Bank 
Rated Capacity = 750 kwh 
Rated Power = 10 kw 

• FLYWHEEL 

Steel Rotor 

Variable Ratio Transmission 
Rated Capacity = 50 kwh 
(24,000 pound weight) 

Rated Power = 20 kw 

• THERMAL 

Parafin Wax Phase Change Media 
Rated Capacity = 500 kwh 
Rated Input Power = 50 kw 

LOGIC 

• Priorities for generated power 

1 - Battery Storage 

2 - Flywheel Storage 

3 - Thermal Storage 

4 - Surplus to Utility 

• Utility Backup for Storage Media 

• Battery Services AC Base Load and 
DC Load, Flywheel Services AC Peak 
Load 
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(a) AC LCAD - NOV. TO APR. (b) AC LOAD - MAY TO OCT. 



(c) DC LOAD (d) THERMAL LOAD 


Figure 3.1-4 Daily Load Profiles for Hybrid Model 
BCS 40259-1 /iff- 49 


profiles used for the electrical and thermal loads. Since the AC load is a 
function of both time of day and month, a two dimension table lookup component 
was used in the system model to represent the load profiles shown in (a) and 
(b). The mean value of the combined electrical, thermal, and control loads is 
10.8kw. It should be noted that this model was chosen simply to illustrate the 
capabilities of the SIMWEST program and is not a proposed concept or economic 
system design. 

In addition to the hybrid model simulations, several variations of this model 
were sized and simulated for comparison purposes. The alternative models use 
only solar photovoltaic generation and wind power generation, respectively. 
Table 3.1-2 shows the system parameters used to size each of the models. 
Generation was sized to the load by increasing rated output until the energy 
surplus exceeded that required from the utility backup to satisfy the thermal 
load. Rated input power to the battery was also adjusted, so that the 
electrical loads are primarily satisfied by solar and wind generation with 
excess power going to thermal storage. The rated solar P-V parameters are 
considerably larger than the Hybrid and Wind model parameters since solar 
energy generation is restricted to an 8-9 hour period each day. 


Table 3.1-2 Test Case Load Matching Parameters 



MODEL PARAMETER 

HYBRID 

SOLAR P-V 

WIND 

• 

! 

Wind Turbine 
Blade Radius 

6.4 m 


8.8 m 

i 

• 

Rated Wind 
Output Power 

20 kw 


36 kw 

• 

Solar Array 
Area* 

189 m 2 

459 m 2 


. 1 

• 

Rated Solar 
Output Power 

30 kw 

75 kw 

— 

• 

Battery Rated 
Input Power 

10 kw 

18 kw 

10 kw 


1 Number of solar cells and collector dimension parameters varied pro- 
portionately. 
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3.2 DAILY SYSTEM PERFORMANCE 


In this section daily and short term performance of the energy storage systems 
is illustrated by analyzing environmental and storage time plots over a weekly 
period. Figures 3.2-1 and 3.2-2 show the collector solar insolation and wind 
speed inputs for the week selected. Both the wind and solar inputs are strong 
during the first three days, and then fall off for several days, gradually 
increasing toward the end of the week. These inputs are typical for this time 
of year. 

The battery, flywheel and thermal energy storage time plots for the Hybrid 
model are shown in Figures 3.2-1 to 3,2-5, respectively. Battery storage 
energy is a rather smooth curve, not responding greatly to the daily fluctua- 
tions in generation and electrical load. The total swing in battery storage is 
in fact less than half of its capacity during this period. Flywheel storage, 
on the other hand, is quite responsive to daily fluctuations in generation and 
swings rapidly to full capacity when generation is strong and to fully dis- 
charged when generation is weak. Thermal storage is not very responsive to 
daily fluctuations but absorbs most of the output when generation is strong. 
Thermal storage tends to become fully discharged during several day periods of 
weak generation. 

For comparison purposes, the energy storage time plots for the Solar P-V model 
are shown in Figures 3.2-6 to 3.2-8. These plots are similar to those for the 
Hybrid model except that the daily fluctuations are more pronounced, and more 
regular than those with wind generation. 

It has been found that performance data using four seasonal week long simula- 
tions is quite suitable for solar systems trade study and design tasks. For 
example, the above analysis of the Hybrid model points out several areas in 
which the system design might be improved: The plots of battery storage tend 
to indicate that the battery is oversized and its capacity could be reduced. 
The flywheel seems well sized to its purpose, but frequently calls on a utility 
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Figure 3,2-3 Battery Energy Storage - Hybrid Model 
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Figure 3.2-5 Thermal Energy Storage - Hybrid Model 
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Figure 3.2-6 Battery Energy Storage - Solar P-V Model 
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Figure 3,2-7 Flywheel Energy Storage - Solar P-V Model 
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Figure 3.2-8 Thermal Energy Storage - Solar P-V Model 



for charging even when the battery has substantial capacity. This suggests 
that the battery back up the flywheel and service its load whenever the fly- 
wheel is fully discharged, rather than calling on a utility. Finally, the 
thermal storage unit is too small to supply a substantially higher percentage 
of the load, since increased generation capacity tends to be surplused to the 
utility. 

3.3 MONTHLY POWER GENERATION 

The yearly wind and solar simulations are best set-up as twelve monthly simula- 
tions in order to generate partial results from printer plots, histograms, and 
numeric tabular output. This section summarizes monthly results from the year 
long simulations. 

Figure 3.3-1 shows the average monthly generation for the Hybrid solar model. 
The wind turbine output is highly variable from one month to the next and does 
not show a significant seasonal trend. This may be an artifact, however, 
resulting from TMY data selection of "typical" meterological months from dif- 
- ferent years of SOLMET (Ref. [3]) weather and insolation data. The solar array 
output, by contrast, is a fairly smooth curve with a distinct mid-summer peak. 
Utility back-up generation is complementary to the wind and solar generation, 
with no back-up required for those months with strong wind generation. 


Figures 3.3-2 and 3.3-3 show the monthly output histograms for wind and solar 
generation during May. This was one of the strongest months for wind gener- 
ation and the wind turbine operated at rated ouput (20kw) about 30% of the 


time. The solar array output, although more consistent than the wind, produces 
substantial output only for a short period each day. Only 20% of the time did 
the solar array produce more than 50% of its rated capacity. Thus the monthly 


plant factor or ratio between average output and rated output was 0.55 for the 
wind turbine and 0.20 for the solar array. 
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Figure 3.3-1 Average Monthly Output Power For Hybrid Model 
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Figure 3.3-2 Histogram of Wind Turbine Output Power - May 
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Figure 3.3-3 


Histogram of Solar Array Output Power - May 
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Figures 3.3-4 and 3.3-5 show the average monthly output for the Wind and Solar 
P-V models. The average output power from the wind turbine in 3.3-4 is higher 
than the combined wind and solar generation for the Hybrid model, but the 
average utility power required is also higher due to the larger power deficits 
during less windy months. The average power surplused is about equal to that 
demanded from the utility in this case. Figure 3.3-5 shows that the Solar P-V 
system has a substantial power deficit during the winter months which changes 
to a small surplus in the summer months. This indicates that the solar array 
may be undersized for yearly system operation. 

3.4 LOAD SERVICE PERFORMANCE 

Many of the component models monitor power flow statistics, such as total 
energy into or out of a component since the beginning of simulation. The use 
of these statistics and histograms enable the SIMWEST user to analyze subsystem 
and component level performance of a wind or solar simulation. This is 
illustrated in this section by analyzing the distribution of power flow to the 
loads, and storage system efficiencies for the Hybrid model. 

Figure 3.4-1 shows the average monthly distribution of input power to storage 
and surplus. The battery input is fairly constant, ranging between 5 and 6kw 
each month except January. Similarly, the flywheel input is fairly constant, 
ranging between 1.5 and 2kw each month. The thermal vessel, and to a lesser 
extent surplus power follows the ups and downs of the wind based generation. 
It should be noted that this figure simply confirms the analysis of daily and 
short term performance obtained with weekly time plots of energy storage, i.e., 
Figures 3.2-3 to 3.2-5. 

The monthly distribution of power flow to the electrical loads is not of much 
interest since the battery supplies the DC load and the flywheel and battery 
supply the AC load in proportion to the peakload and baseload demand, 
respectively, except for a small amount of power supplied by a utility when the 
flywheel is discharged. Figure 3.4-2 shows the monthly percentage of power 
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Figure 3.3-4 Average Monthly Output Power- Wind Model 
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Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
Figure 3.3-5 Average Monthly Output Power - Solar P-V 
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Figure 3.4-1 Average Monthly Power Distribution 
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Figure 3.4-2 Monthly Thermal Load Service - Hybrid Model 
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to the thermal load from solar and wind generation, the percentage from a 
utility, and the percentage of tfiermal load not met. Solar and wind generation 
supply at least 75% of the thermal load each month except July and September 
when wind generation is low. This suggests increasing the solar generation 
capacity, relative to the wind in order to reduce the utility supplied load 
during the summer months. The small percentage load not met is due to an 
assumed temperature deadband when thermal storage temperature drops below a 
minimum useful level. 

Table 3.4-1 summarizes the storage system performance in meeting the loads. 
The values shown represent total subsystem level performance, i.e., input 
power to the flywheel is the DC power to the motor and output power is the power 
delivered to the loads by the generator and inverter. The battery efficiency 
is higher than would be expected for a conventional lead-acid battery, but this 
is easily corrected by appropriate scaling of the battery terminal resistance 
parameter. This would lower the overall storage efficiency and thus require 
greater total generation to meet the load. 


Table 3.4-1 Average Storage System Statistics 



Input. 

Power 1 

Output 

Power 

Turnaround 

Efficiency 

Battery 

5.6 

5.2 

0.94 

Flywheel 

1.8 

1.3 

0.75 

Thermal 

5.9 

4.2 

0.71 

Total 

13.3 

10.8 

o.8i ■■■■;■ 


^ Adjusted for net change in energy storage 
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3.5 SUMMARY 


The Hybrid solar test case simulations have illustrated the general ability of 
the SIMWEST program to model and analyze complex solar and wind system models. 
These studies have shown that using four seasonal week long simulations rather 
than expensive year long simulations gives an adequate data base for system 
level design and trade study analyses. Year long simulations are then 
recommended for final confirmation of the selected system designs. For 
simplicity, we have neglected cost economics constraints associated with an 
actual system application. However, the economic and life cycle cost analysis 
described in Section 1 could have been used to investigate system economics. 
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4.0 TECHNICAL APPROACH FOR POWER DYNAMICS LIBRARY 


4.1 BACKGROUND 

The technical approaches studied for development of a dynamic component li- 
brary to model transient effects on generation, distribution, and storage of 
power are summarized in this section. The objective of this library develop- 
ment was to provide component models for analysis of transients such as wind 
gusts, transmission line faults and generation or storage switching on an 
electrical network. Only those phenomena with physical time constants on the 
order of hundredths of a second to ten seconds were envisioned as within the 
scope of the component modeling. This is the primary region of interest for 
power system transient stability analyses. Most of the dominant machine time 
constants for power generation and load distribution fall within this region. 
The component models formulated do not model detailed control and design phe- 
nomena such as wind turbine structural vibrations or synchronous generator 
excitation voltage modeling. Similarly, system phenomena with time constants 
on the order of minutes to hours were not modeled, except to the extent 
necessary to specify steady state operating conditions. 


One of the primary considerations in formulating the dynamic library was to 
make the codes easy to use. Simplified dynamic models were envisioned for this 
purpose so that a user could create and debug a simulation quickly, and run 
simulations with reasonable cost expectations. Other user related factors 
such as model generality, documentation clarity, output capability and model 
adaptability were also considerations during the model definition effort. 
Within this framework, a primary objective of the model definition effort was 
to obtain component models whidh could be easily modified or augmented to 
incorporate control system interaction with the network, i.e., excitation reg- 
ulation on a synchronous generator. A user would build up a control system 
from transfer function, algebraic operations, saturation* and other library 
components and link the control system to the affected physical components. 
Such control system modeling would be desirable, for example, to provide 
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damping for electrical machine transients or lightly damped DC transmission 
lines with series inductors. In addition, control system modeling is important 
for evaluating system stability under fault current and over voltage con- 
ditions. 

The analysis capabilities of the EASY program for control system design and 
evaluation were to be incorporated into SIMWEST as part of the dynamic library 
development effort. (Table 4.1-1 summarizes these analysis capabilities.) The 
proposed marriage of the EASY program analysis capabilities and components for 
power generation, transmission, and distribution would provide a unique re- 
search tool for dynamic analysis of solar energy/util ity/energy storage sys- 
tems. The resulting program would thus provide the capability of evaluating 
control system designs as well as performing transient stability analyses. 

During the formulation of the dynamic component library it became evident that 
initializing the user's model to steady state power flow conditions is a more 
difficult technical problem than originally anticipated. In general, when 
simulating a dynamic system for the purpose of investigating stability, steady 
state operating conditions are a prerequisite for investigating system tran- 
sient behavior. The EASY program has a steady state analysis capability which 
yields a steady state configuration given user specified initial conditions. 
However, this capability is not sufficient for steady state power flows since 
there are voltage and power angle operating constraints which are necessary 
when modeling power system operation. For a system with more than a modest 
number of busses, or for a user not familiar with power system behavior, 
specification of appropriate initial conditions to balance real and reactive 
power within these constraints could be an arduous task. Consequently, a 
technical approach for initializing the user's model was formulated which 
employs a classical load flow analysis, and a separate set of parameter compu- 
tations at each network bus. An EASY steady state analysis would also be 
required in order to initialize the dynamics of other components such as wind 
turbines and utility steam driven turbines. 
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TABLE 4.1-1 EASY PROGRAM ANALYSIS CAPABILITY 

ANALYSIS PURPOSE 

• SIMULATION CONCEPT VERIFICATION. SYSTEM OPTIMIZATION AND DESIGN 

• MODEL LINEARIZATION _JL.JL_ — — 

• STEADY STATE CALCULATIONS INITIALIZE TRANSIENT SIMULATIONS, INVESTIGATE STEADY STATE 

CONDITIONS AND SYSTEM STABILITY 

• ROOT LOCUS I 

• STABILITY MARGIN > CONTROL SYSTEM DESIGN, SYSTEM STABILIZATION AND SYNTHESIS 

• OPTIMAL CONTROLLER I 

• EIGENVALUES I 

> FREQUENCY DOMAIN ANALYSES, SYSTEM STABILITY 

• TRANSFER FUNCTION I 


-vj 

CO 



This technical approach for initializing a transient power simulation was 
considerably beyond that originally anticipated. Since the estimated increase 
in funding to implement this approach exceeded the project scope, further 
development in this area was terminated. It was decided, however, to document 
the technical approach and component models formulated, since the capabilities 
envisioned through development of such a program would provide a unique re- 
search tool for investigation and design of innovative power system generation 
and storage technologies. 

The SIMWEST dynamic component models which were studied in depth during the 
model formulation effort included a load flow component, network solver, syn- 
chronous machine models, induction machine models, and a wind turbine compo- 
nent. Descriptions of these components and the general problem formulation 
summarized above are provided in subsequent sections. These components de- 
scribe all the network related components originally contemplated, except for 
an AC/DC converter model. The simplified converter model described in refer- 
ence [5] appears to be compatible with this problem formulation. However, the 
model formulation effort was terminated before this model could be evaluated 
for inclusion in the dynamic library. 

4.2 GENERAL PROBLEM FORMULATION 

A power system model for transient stability simulations consists of the physi- 
cal machines or subsystems which generate power into an electrical network, the 
transmission grid which allocates power flow between the network busses, and 
the machines or subsystems which constitute the loads. From a modeling view- 
point, the power network consists of a set of busses or nodes which are 
connected by transmission lines, and which is characterized by (complex) volt- 
age and power at each node. A power simulation consists of four basic calcula- 
tion steps: 

1) Calculation of the machine dynamics or differential equations, given the 
current state of the power system. 
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2) Numerical integration of the system dynamics to the next time step. 

3) Calculation of two of the nodal electrical states for each bus, based on 
the current machine dynamics. 

4) Solution of the electrical network equations for the remaining two elec- 
trical states at each bus. 

Steps 3) and 4) may involve iteration, depending on the complexity of the 
system models. The component models which were selected for the dynamic 
library do not require iteration and thus would economize on computer time for 
transient simulations. 

Three types of busses are modeled in the network component which computes step 
4) above: 

• Busses which provide a specified current injection into the network. 
This includes machine busses such as generators, motors and con- 
verters . 

• Reference or infinite busses which specify a known bus voltage and act 
as infinite power sources or sinks and thus stabilize the network. 

• Internal or passive busses used to provide power flow branching in the 
transmission network and to represent loads. 

Figure 4.2-1 shows a power network model illustrating all of the above types of 
network busses. The reference bus is used to model an interconnection with an 
external utility network. The branching bus is used to furnish redundant 
transmission capability to the loads. The other busses all inject current into 
the network and result in positive power flow for generation and negative power 
flow for the loads. 
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REFERENCE BUS 


Figure 4.2-1 Example Power Network Schematic 





The general flow of power calculations for the components formulated in this 
section is illustrated by Figure 4.2-2 for the five bus model of Figure 4.2-1. 
Following integration of the system dynamics, the first calculations in the 
model are the Norton current injections into the network (blocks (l) and (?)) 
for the induction and synchronous machine busses. The reference bus furnishes 
a constant voltage into the network component. The internal bus and the two 
load busses have been replaced by an equivalent circuit in the network and thus 
are transparent to the solution process. The network component (block (3)) 
solves the circuit equations for bus voltage at all current specified busses 
and for current flow at the reference bus. The wind turbine-generator 
components (blocks ©, © and (6) ) then compute the mechanical torque into 
the generator and the differential equations of the machine dynamics. (Feed- 
back connections are denoted by dashed lines in Figure 4.2-2.) Finally power 
flow and machine dynamics of the motor-transmission-flywheel components are 
computed (blocks © and ®). 

At the beginning of each simulation the system model is initialized by a load 
flow computation for specified steady state power flow at each bus node, and by 
a steady state analysis to obtain initial conditions for the machine dynamics. 
The load flow output and initialized state variables are stored for use in 
later simulations or later runs. During a given simulation run the network 
topology remains constant, i.e., no switching in or out of generators, loads, 
or transmission lines. The user introduces such transients by running a 
simulation prior to the switch, then beginning a new simulation with changes in 
the network, and running a subsequent simulation to analyze power transients. 

4.3 LOAD FLOW INITIALIZATION 

A load flow calculation is a nonlinear equation solution method for obtaining 
the steady state real and reactive power at each nodal bus. This calculation 
procedure is widely used in the electric utility industry, and is the standard 
method for initializing transient stability programs. Its purpose in the 
problem formulation above is to provide appropriate shunt constants for the 
load busses and to provide initial conditions for synchronous and induction 
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Figure 4.2-2 Example SIMWEST Model Schematic 




machine injections. If one were to use the EASY program steady state analysis 
and nominal shunt parameters to balance real power flow, it might be possible 
to obtain steady state power flows without a load flow calculation. However, 
such solutions would generally have poorly balanced reactive power generation 
resulting in large voltage magnitude differences throughout the network. Thus 
a load flow calculation is essential in order to balance both real and reactive 
generation to the loads. 

Three types of busses are modeled in load flow calculations: 

1) load busses (including induction machines) where the net real and re- 
active power injections P and Q are specified and the load flow solves for 
voltage magnitude |V ( and angle 5, 

2) Voltage controlled busses such as synchronous machines where the voltage 
magnitude Jv j and real power P are specified and the load flow solves for 
reactive power Q and voltage angle 8 , 

3) Reference busses where Jv and 8 are specified and the load flow computes 
power injections P and Q which balance real and reactive power flow 
throughout the network. 

Passive busses may be ignored since a network reduction may be performed to 
eliminate such busses prior to the load flow. 


Mathematically, the load flow is a procedure for solving the n nodal power flow 
equations 

: s ,n - ^ 

where the (complex) currents I m satisfy the network equations 



“ 2 k Y mk 



(4.2) 
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and (Y mk ) denotes the complex valued admittance matrix of the transmission 
network. Using the polar representations 



- u 4 A 
* 'mk 



e J0mk^ 


and substituting (4.2) into (4.1) yields 


S m = \ V m V k Y mk| e yW + 6 k 8 m^ (4.3) 

Separating (4.3) into real and imaginary parts, one obtains 

°-.V S k |v m V k V mk j cos <e mk + 6 k - y (4.4) 

0 = " 2 k v m v k V mk| s1n (e mk + 5 k ' (4-5) 

Equations (4.4) and (4.5) represent a nonlinear system of equations of the 

general form 


0 = g(x,y) 


(4.6) 


where x represents the 2n unknowns to be solved, y represents the 2n nodal 
conditions which are specified, and g represents the 2n equations (4.4) and 
(4.5). 


Given an initial estimate x Q , (4.6) may be solved by Newton-Raphson iterations 

-|f (x is y) • 6x i = g(x i ,y) ;v ' (4.7) 

Xj +1 - x. - 6x. (4.8) 

Although a number of advanced techniques exist for efficiently updating the 
Jacobian matrix and for solving the linear equations (4.7), these methods 
are probably not warranted for the small dimension networks envisioned for 
SIMWEST transient stability simulations. Consequently, the Jacobian matrix 
Would be evaluated each iteration with the particular structure of equations 
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(4.4) and (4.5) used to reduce the evaluation effort, and Gaussian elimination 
then used to compute . 

Once a load flow solution to (4.4) and (4.5) is obtained, appropriate network 
constants and bus states may be supplied to the network and machine components. 
The current injections for each bus may be obtained by solving (4.1), i.e., 

V; s V v ; . - (4.9) 

Similarly, the shunt impedance for constant load busses satisfies V = Z s I, and 
thus 

S = V*I* = 1 1 1 2 * Z s 

or equivalently A 

z s « <P + JQ)/| (4.10) 

After these quantities are computed, the Y admittance matrix, the nodal current 
and voltage quantities I ,V m and the shunt impedances Z s are output to disc 
storage for later use in initializing transient simulations. 


In order to fully initialize a transient simulation, the initial conditions of 
all state variables must be selected consistent with the steady state load 
flow. This is easily done using the EASY steady state analysis, but requires 
adding new state variables to the model whenever a specified load flow condi- 
tion; is the output of other components. For example, in the model illustrated 
in Figure 4.2-2, the output power of the generator is a function of wind speed 
whenever V UTWn is less than rated wind speed. In this case a feedback control- 


ler is used during the steady state analysis to obtain the value of 
corresponding to the specified output power P GEN (See Figure 4.3-1). The 
steady state routine will iterate until P G converges to P G£ r N , solving for the 
corresponding v W j ND value. This part of the model can be switched off after 
model initialization so that is independent of the output power P G . 
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Figure 4.3-1 Control Schematic to Initialize Example Model 
4.4 NETWORK COMPONENT MODEL 

This component provides solutions for the voltage and power flow at each bus, 
given the injected current for machine busses and the reference voltage for 
reference or infinite busses. On the first pass through the model during a 
simulation or steady state analysis, the network admittance matrix is con- 
structed and a triangular decomposition is performed. The admittance matrix 
may be consructed either from the load flow matrix, or from user specified 
input data, and includes a shunt impedance for load busses and machine busses. 
During each subsequent pass through the model the linear network equations are 
solved and bus power and voltage variables are output. 

1) Admittance Matrix Calculations (First Pass Only) 

The network admittance matrix may be constructed either from the load flow 
admittance matrix or from user specified inputs: 


• Load Flow Initialization 


Read in the load flow admittance matrix Y and shunt impedances Z m for all load 
busses. Recompute the diagonal elements of Y such that 


Y* 

mm 




82 


BCS 40259-1 



User Specified Initialization 


V lm = W- 

Y ™ ■ - E l*n Y l m + (>C, m /2) + 1/Z m 

where 

Z lm = Impedance of the transmission line between nodes 1 and m 

w C lm = Line capacitance to ground between nodes 1 and m, 

multiplied by synchronous frequency (120^ 
radians/second) . 

Z m = Shunt impedance of a load or machine bus. 


2) Network Reduction (to eliminate passive busses) 

Partition the admittance matrix into external or injection busses and 
internal or passive busses: 

. ' 1 " ' ' . • 

•. ■ . ’ I • ‘ 



The equivalent reduced network for the external busses is given by 



y ei Y ii 




Where the matrix X=Y E jYjJ is obtained by Gaussian elimination of the linear 
matrix equation 



3) Solution of the Network Equations 

Partition the admittance matrix Y' into current injection busses denoted 
by the symbol P and reference busses denoted by R: 
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Let Ip, I R and V p , V R denote the complex valued current and voltage 
vectors corresponding to the current injection busses and the reference 
busses. Then the network equations may be written as 



Y pp Y pR 


yT y 
y PR y RR 



(4.11) 


We may expand (4.11) into partitioned equations for I D and I D and formally 
solve these equations for V p and I R : 

Vp = VpJ (ip - y pR V R ) (4.12) 

*R ■ V PR '*'? + V RR V R (4.13) 


It is more efficient and ^accurate to solve (4.12) using a triangular 
factorization of Y p p : 

Ypp = L * U (4.14) 

This factorization is performed during the first pass and then stored to 
solve (4.12) by forward and backward substitution, i.e. , we solve the two 
sets of equations 

LX = Ip-Yp R V R (4.15) 

U\T p = X (4.16) 

Once Vp is obtained, is computed from (4.13). The complex power flow at 
each bus is then computed from the bus current and voltage variables: 
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(4.17) 


P m + jQ m = V m I* 
m °^m mm 

4.5 SYNCHRONOUS MACHINE MODEL 

The synchronous machine model selected for the SIMWEST dynamic library is a 
compromise between the simple classical machine which models only the machine 
rotor mechanics, and more detailed machine models which may include excitation 
control, rotor transient and subtransient dynamics, and stator transients. 
Specifically, rotor electric and mechanical dynamics are modeled, field volt- 
age is constant, and stator and subtransient dynamics are ignored. The model 
is sufficiently detailed to enable a user to add excitation or governor con- 
trols to the system model, and is simple enough to minimize the computations 
required to solve the network equations. The synchronous model may be used to 
represent either a motor or a generator, depending on the algebraic sign of 
injected power into the network. Another component is also needed with the 
synchronous machines to compute the Norton equivalent current into the net- 
work. These models are described below. 

4.5.1 Synchronous Machine Model 

The model given here is a simplified version of the machine model in reference 
[6]. The model is most easily described in terms of the machine direct and 
quadrature rotating coordinate system (d-q axes). Four state variables are 
used to represent the machine dynamics: 

5 = rotor angle in radians 

w = rotor shaft speed in radians/second 

e'd» e'q = d-q axis representation of the machine internal 

voltage E 1 (See Figure 4.5-1) 
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Figure 4.5-1 Network and d-q Axis Coordinates 
The basic equations for this model are: 

1) Stator Electrical Equations 



where 


*d’ : q 

Sr 

bus injected current 

in d-q axes 

V d> \ 

. ; = 

bus terminal voltage 

in d-q axes 

«s 

= 

stator resistance 


x 'd- 

= 

transient reactance 


Machine Dynami 

cs 


t e 

= 

< e 'd • : d * °'d • *«>* 

• ■ 

Jw 

= 

t m " t e " D ' 



6 = a) - w 

0 


where 

T e = electrical torque output 


(4.18) 


(4.19) 


(4.20) 
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0 = rotor drive shaft inertia 

D = machine damping 

Tjvj = mechanical input torque 

u> 0 = synchronous frequency (120 radians/second) 


3) Internal Voltage Dynamics 




where 



UqU,, - x' d ) - e' d )/T' q 

( E f - e ', - 'd< x d * X 'd» /T 'd 

= synchronous reactance 

= open circuit transient time constants 

= field voltage of excitation system 


(4.21) 


During a simulation the terminal voltage V from the network is transformed to 
d-q axes using the rotation equations: 


V d = sins • Re (V) - cos6 • Im (V) 

Vq = cos <5 • Re (V) + sin<5 • Im (V) 


(4.22) 


Equations (4.18) and (4.19) for current and electrical torque are then solved 
and the differential equations (4.20) and 4.21) are computed. 


• Initialization (First Pass Only) 

When load flow initialization is desired the state variables w and 6 are 
computed using 


(D 

tan 5 


(Im (V) + R $ • Im (I) + x d • Re (I))/ 
(Re (V) + R $ * Re (I) - x d *Im (I)) 
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where I, V denote the bus current and terminal voltage from the load flow. 
Given the angle , the d-q axis components for I and V may be computed using the 
rotation equations (4.22). Initial values of e* d and e' q are then obtained 
from (4.18), i.e., 

e 'd ■ V d + V ‘d-xV ‘q 


The field voltage and mechanical torque T M values driving (4.20) and (4.21) 
to zero are then computed: 

T M = !e 'd • r d + e 'q • V /o> o + 

E f = e 'q + V < x d ' x 'd> 

If the machine is being used as a generator, the input torque may require 
initialization by the procedure indicated in the last part of Section 4.3. On 
the other hand, if the model is being used as a synchronous motor then may be 
a state variable feedback from a downstream component. In this case the above 
computation for T^ initializes this state. 

The admittance matrix Y* for the network is also modified to represent the 
Norton equivalent stator circuit, i.e., 

Y' = Y mm + 1/(R_ + jxi) 
mm mm s ° a' 

where m is the index of the synchronous machine bus. 

4.5.2 Current Injection Component 

The current injection component computes the Norton current I into the network 
given the state variables E' and S. The stator equations (4.18) neglecting V d 
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and Vq are used to solve for d-q axis components I d and I q and then a trans- 
formation is performed to obtain the synchronous frame components. Specifi- 
cally, solving (4.18) for I d and I q yields 

I d = < R S • e d + x 'd B 'q> /(R s + x 'd • x 'd } 

■q ■ <- x 'd e 'd + R s e 'q> /(R s + x 'd • x ‘d> 

and thus 

T re = +s1n 6 * 1 d + cos 6 * r q 

= - coss - r d + sin 5 

4.6 INDUCTION MACHINE MODEL 


The induction machine model is a simplified model which neglects stator tran- 
sients and includes' rotor electrical and mechanical dynamics. The machine 
model may be used to represent either a motor or a generator depending on the 
algebraic sign of the input torque. A current injection component is also 
required to compute the Norton equivalent current into the network. 

Three state variables are used to represent the machine dynamics: 

to = rotor shaft speed in radians/second 

0 I p • 

d’ q = d-q axis representation of the machine internalvolt- 

age E’. . : 

The axes are chosen such that the direct axis coincides with the network frame 
real axis and the quadrature axis coincides with the network frame imaginary 
axis. Torque and slip are positive for a motor, negative for a generator. 

The basic equations for this model in d-q axes are 
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1) Stator Equations 


x' - x s - x M /x r 


e 'd- v d\ hs *• \/ i d 


e YV V*' - R s/Vq 


(4.23) 


where 


v y 

y i,*. 


stator reactance 

magnetizing reactance between stator and rotor 
rotor reactance 
stator resistance 

d-q axis components of bus terminal voltage V = V q 

d-q axis components of injected current I = I d +jlq 


2) Machine Dynamics 


< e d ‘ J d + e 'q ‘ V/“ 
T e - T m - Du 


(4.24) 

(4.25) 


where 


electrical torque 
rotor drive shaft inertia 
machine damping 

mechanical torque (input variable) 
synchronous frequency (120 it radians/second) 


3) Electrical Dynamics 


K - W) e 'q - "f (e 'd + K - *') X q> 


(4.26) 
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■ (“ - %) «'d - -r < e 'q - ( x s- x ') I d > (4 - 27) 

where 

x r 

T 0 = , R f - rotor resistance 

During a simulation, (4.23) is solved for I. and I , and then equations (4.24) 
to (4.27) are used to compute the differential equations for w, e'^ and e'^, 

• Load Flow Initialization (First Pass Only) 

Given the load flow P +jQ and bus voltage V, the steady state slip or rotor 

speed corresponding to real power output is computed. Then a shunt reactance 
balancing reactive power flow is obtained. Finally, I d and I q are obtained in 
order to compute the initial states e' d and e' q . The equations below are based 
on the circuit model shown in Figure 4.6-1. 



Figure 4.6-1 Steady State Circuit for Induction Machine 
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1) Steady State Slip 


C 2 s 2 + C x s + C Q = 0 (4.28) 

where 



The classical quadratic formula may be used to solve (4.28) for s. w is then 
given by 

(1 - s)io 0 \ 7/;: (4.29) 

2) Shunt Reactance and Admittance Matrix Update 

1 Q e = ^2- ( x s ( R r /s2 + 4 ) * V x m) ; 7 < 4 - 30 > 

with 

V> j R s V s • x s x r + X M ; 

B * R $ x r + R r x $ /s 

P e +jQ e represents the power flow from the induction machine. Since the load 
flow power Q is in general not equal to Q e , a shunt reactance is used to balance 
reactive power : 

j(Q - Q e ) = v • (v*/ - jx Q ) = j lv| 2 /x 0 
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Thus 


X 0 = 1 v| 2 /(Q - Q e ) (4.31) 

The admittance matrix Y for the network is then modified such that 

Y' = Ymm + 1/{R„ + JX ' ) - j/X„ 
mm mm ' v s ' o 

3) Initial Voltage States and Torque 

The current components I d and 1^ are obtained by solving 


(P e + JQ e ) = V* I* 



e' d and e' q are then obtained from (4.23). Mechanical Torque is given by 

t m = t e - Dai 

where T^ is computed from (4.24) 

4.6.1 Current Injection for Induction Machine = . ::r , 

The current injection to the network for the Norton equivalent of the induction 
machine is obtained by solving (4.23) for I d and I , neglecting the bus voltage 
terms : ' 

: ; I d = -(V e' d + x 1 e‘ q )/{R| * (x') 2 ) 

Iq ’ (x* • e' d - R s - e- q )/(R| + (x-) 2 ) 
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4.7 WIND TURBINE MODEL 


The wind turbine model is a simplified model of the turbine dynamics and blade 
pitch controls, ignoring the structural dynamics of the tower and flexing 
motions of the rotor blades. This model is based on a wind turbine simulation 
developed for studying pitch control effectiveness of the Mod-2 wind turbine 
under development for DOE. The model includes a twice-per-revolution 
disturbance due to passage of a blade through the tower wind shadow, and 
damping of the rotor hub oscillations by use of a flexible shaft or hydraulic 
coupling to the generator. 

Figure 4.7-1 shows the main components and variables for the wind 
turbine/generator system. The blade aerodynamics are summarized by 
coefficient of power (C ) curves as functions of wind tip speed ratio (*) and 

r 

blade pitch (g). The torque output at the hub (T^yg) is computed based on C p 
and a blade angle dependent torque disturbance. The shaft dynamics consist of 
differential equations for rotor reference angle and speed at the hub (0y, ^), 
differential equations for shaft flexing angle and speed ( 5 S » W S ) at the gearbox 
coupling, and calculation of output torque (T g ) to the gearbox. The gearbox 
model is a simple table lookup used to compute power loss in the output torque 
(T m ) to the alternator. The pitch controller is a linear control system for g 
which includes integral and proportional feedback of hub speed error and 
generator output power error, and a notch filter to attenuate twice-per- 
revolution disturbances. Basic equations for these components are outlined 
below:. ' 

• Blade Aerodynamics 

Given and the tip speed ratio A is computed: 

= r * W H /V WIND (4.34) 

where r = rotor blade radius, 
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Rotor Speed 



Model-Schematic 





The wind power generated at the hub is then given by 


P WIND = c p( e » x ) ,7rr * pV wiND //2 


where p = air density 


C p (3,x) is obtained by interpolation of C p curves versus X for fixed B. 
(See Figure 4,7-2). 



Figure 4.7-2 Coefficient of Power Design Curves 

The torque output at the hub is the sum of the torque obtained from (4.35) and a 
blade angle dependent disturbance: 

T HUB " + s i n ( 2 e H + T ) P WIND^^H (4.36) 

where. 


K 

e H 


disturbance gain (nominal value =0.1) 
rotor hub angle relative to horizontal blades 
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• Shaft Dynamics 

The shaft dynamics are represented schematically by Figure 4.7-3. The 
differential equations for <$ s and w s are: 

0 ^ = ^ (4.37) 

J H W H = T HUB ' °H W H + D Q^s"“H^ + K Q 5 s ^.38) 

6 $ = o) s - - (4.39) 

J s “s = " D s “s + D q( w H “ “s^ ' ,C Q 6 s + D 6 ^“g ■ “s^ (4.40) 

where 

Jh»J s = rotor hub and drive shaft inertia 
D h ,D s = hub and shaft damping 
Dq,Kq = quill shaft damping and compliance 
0 G = hydraulic damping between shaft and generator 
= generator speed/gear ratio. 


The output torque to the gearbox is given by 



Figure 4.7-3 Turbine Shaft Dynamics 
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• Gearbox 


The power loss in the gearbox is modeled by a lookup table L (T) : 

T m = G-T s - L(G-T S ) (4.42) 

where G denotes the gear ratio. 

• Pitch Controls 

Representative pitch controls are easily summarized using transfer function 
equations. A raw pitch command 3 C is generated using a fixed value or integral 
plus proportional feedback: 

3 C - (Kj + K 2 /s) (P REp - P 8 ) (4.43) 

+ (K 2 K 3 /S) ■* 

where P RFF is the rated output power of the turbine/generator and w REF is the 
reference speed for synchronous operation. The pitch command is notch filtered 
and input to a pitch actuator. The pitch output 3 is given by 

Actuator Notch Filter 

3 = (l/(T A s+l)) • (s 2 + .02s w RE p + w r EF )/(s 2 + 0.2s “> REF + w 2 £F ) (4,44 ^ 

where T« = actuator time constant. 

The above pitch controls would not be built into the wind turbine component 
since they are easily constructed from transfer function components, and more 
general ity is obtained this way. 

4.8 SUMMARY CONCLUSIONS 

A general technical approach for development of a transient power system li- 
brary using the SIMWEST/EASY computer program has been outlined. Although the 
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problem formulation is an adaptation of power industry transient stability 
program models and solution methods, its purpose and capabilities are rather 
different. The SIMWEST dynamic component library would be an excellent re^ 
search tool for investigation and design of innovative or unconventional en- 
ergy systems interconnected to a power grid, whereas transient stability pro- 
grams are usually production codes used to represent existing or proposed power 
systems networks. As a consequence of the design objectives SIMWEST has more 
emphasis on user simplicity, modeling flexibility, and control system analysis 
than does a transient stability program. It thus would provide a unique 
capability for solar energy and energy storage power systems research. 


One; conclusion that became clear in the formulation of the dynamic library is 
that an expanded version of the EASY program (version 5) should be used as the 
basis for a power transient program, rather than version 4, which was the basis 
for SIMWEST, Version 5 has vector and matrix input/output and module connec- 

i I 

tion capabilities not available with the earlier program versions. These 
capabilities enable a user to model medium size networks rather than being 
limited to ten active busses, provide more user input and output capability for 
complex models, and simplify development of library component subroutines. 


A major technical problem that arises in creation of transient power system 
codes is that the system models are in general implicit, and require iteration 
techniques to obtain convergence of all system variables. The technical ap- 
proach used has been to formulate the models so that the network solution is 
explicit, reducing the compexity and expense of the solution process. In the 
initialization process, however, a load flow component is needed to solve the 
implicit power flow conditions which balance real and reactive generation to 
the loads. This approach is easily modified for more complex system models 
with implicit variables and nonlinearities by use of the iteration method 
currently used in SIMWEST. Thus, the dynamic library models formulated above 
would provide an excellent base of power system models for transient analyses 
which can accomodate growth and expansion for general solar energy research 
studies. 
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